Prom the Hosftadter to the Fibonacci butterfly 
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We show that the electronic spectrum of a tigth-binding hamiltonian defined in a quasiperiodic 
chain with an on-site potential given by a Fibonacci sequence, can be obtained as a superposition of 
Harper potentials. The electronic spectrum of the Harper equation is a fractal set, known as Hos- 
fateedter butterfly. Here we show that is possible to construct a similar butterfly for the Fibonacci 
potential just by adding harmonics to the Harper potential. As a result, the equations in reciprocal 
space for the Fibonacci case have the form of a chain with a long range interaction between Fourier 
components. Then we explore the transformation between both spectra, and specially the origin of 
energy gaps due to the analytical calculation of the components in reciprocal space of the potentials. 
We also calculate some localization properties by finding the correlator of each potential. 

PACS numbers: 71.23.Ft, 71.23-k, 71.70.Di, 71.23 An, 73.43.-f, 73.43.Nq 



I. INTRODUCTION 

Although the discovery of quasicrystals [1], which are 
alloys with an structure neither periodic, nor disor- 
dered, launched an extensive investigation on quasiperi- 
odic Hamiltonians, in fact the study of such Hamiltonians 
goes back to the old Frenkel-Kontorova model [2] and to 
the pioneer research made in the seventies [3, 4] . Here the 
word quasiperiodic means that in the system there are in- 
commensurate periods; and as a result, the dimension of 
the Fourier space is always bigger than the dimensional- 
ity of the system. One of the most famous quasiperiodic 
Hamiltonians was obtained by Harper in connection with 
a problem proposed by Peierls [3]. The idea was to find 
the spectrum and the wave-functions of an electron in a 
square lattice with a perpendicular magnetic field. Two 
periods are involved in the problem, the electron motion 
in the lattice and the cyclotron frequency [3]. The spec- 
trum as a function of ratio between these periods turned 
out to be a complex set known as the Hostacdter but- 
terfly [3] (see figure 1). Since then, the Harper model 
has been very useful to investigate the transition from 
localized to extended eigenstates, as the spectrum pass 
from pure point to continuous [3] [4] [5]. Between both 
limits, there is a new type of spectrum which is known 
as singular continuous [4] . The corresponding eigenstates 
are called critical and show self-similar properties. For 
certain parameters of the Harper equation, the distribu- 
tion of level spacings follows an inverse power law [6], 
which is a new type of spectral statistics [6], explained 
as a level clustering tendency [7] . ft has been possible to 
find analytical expressions for the wave-functions using 
quantum groups [8]. More recently, the quantum phase 
diagrams [9] and the electronic correlation effects have 
been analyzed [10]. 

Another quasiperiodic system that has been exten- 
sively studied is the Fibonacci chain (FC). This chain 
is the simplest model of a quasicrystal [11]. The impor- 
tance of this Hamiltonian arises because the nature of the 
physical properties of quasicrystals is still not well under- 



stood [12] [13] [14]. Even in the theoretical side there is 
a lack of understanding in how electrons propagate, spe- 
cially in two and three dimensions [15]. As is well known, 
a periodic potential satisfies the Bloch's theorem, which 
tells that the eigenstates of the Schroedinger equation are 
plane waves of delocalized nature, and the energy spec- 
trum is continuous [16]. For disordered systems, like in 
the one dimensional (ID) Anderson model, all the states 
are localized corresponding to isolated eigenvalues [17]. 
In more dimensions, there is a mobility edge which sep- 
arates extended from localized states [17]. For most of 
the quasiperiodic systems in ID, the spectrum is neither 
continuous nor singular, instead a new type of spectrum, 
called singular continuous is obtained [18]. This kind of 
spectrum is similar to a Cantor set, and presents a mul- 
tifractal nature. The corresponding eigenfunctions are 
called critical, and also show self-similarity and multi- 
fractality. In two and three dimensions, the nature of 
the spectrum is not known, although there seems to be 
a kind of mobility edge [19] [20]. However, even in ID, 
where large amount of work has been done, there are 
many unsolved questions, like the nature of conductivity 
[21] or diffusivity [5], the spectral statistics and the shape 
of many of the eigenfunctions [22] [23]. Even in the FC, 
there are no analytical expression for the wave functions, 
except for few energies [11]. 

The Harper and Fibonacci potentials share many char- 
acteristics; for example, both can present a multifractal 
spectrum with self-similar wave functions, although the 
Harper equation can also present pure point and con- 
tinuous spectrum. An interesting question is why Fi- 
bonacci does not show pure point or continuous spec- 
trum. Furthermore, is it possible to build the analogous 
of the Hostadter butterfly for the FC? An understanding 
of these questions will serve to give a better picture of the 
electronic properties of quasiperiodic systems. For exam- 
ple, it can suggest a way to construct analytical solutions 
for the FC in terms of those solutions from Harper. 

In this article, we show that in fact, the Fibonacci 
potential can be made from a superposition of Harper 
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potentials. Then, we can follow the transition of the Hof- 
stadter butterfly to the equivalent in the Fibonacci case. 
This allows to explore the equations in reciprocal space 
of the FC. The layout of this work is the following, in 
section II we show how to obtain the Fibonacci potential 
in terms of Harper. Section III is devoted to a discussion 
of the corresponding spectra using the properties in re- 
ciprocal space, while section IV is devoted to study the 
localization in terms of the correlators of the potential. 
Finally, in section V the conclusions are given. 



II. THE FIBONACCI AND HARPER MODELS 

As a general model we will use a tight-binding Hamil- 
tionian of the type, 

(E-V(n))ll) n = t n tpn+l+tn-lll>n-l, (1) 

where ip n is the wave-function at site n, t n is the reso- 
nance integral between sites n and n+ 1. For the present 
purposes, t n is set to 1 for all sites. V(n) is the atomic 
on-site potential and E are the allowed energies. The 
generalized Harper equation is obtained when [3], 

V(n) = V H (n) = 2Acos(27r^n + u), (2) 

where A > is the strength of the potential, (f> is a pa- 
rameter that contains the ratio between the electron cy- 
clotron frequency and the elementary quantum flux, and 
v is a phase shift. For a rational <fr, Eq. (1) can be 
solved by Bloch's theorem since the problem is periodic, 
although the value of such theorem is very limited since 
the coefficients in Fourier space of the solution form a 
dense set [3]. For </> irrational, the spectrum depends on 
the value of A. For A < 1, the spectrum is continuous 
with extended wave functions, for A > 1 the spectrum is 
made from pure points and localized solutions. For A = 1 
the spectrum is singular continuous with self-similar wave 
functions. 

The other potential that we will consider, is the sim- 
plest model of a quasicrystal, obtained when V(n) has 
two possible values, Va and Vb following the Fibonacci 
sequence (FS). The FS is build as follows: consider two 
letters, A and B, and the substitution rules, A — > B, 
and B — > AB. If one defines the first generation se- 
quence as J 7 ! = A and the second one as = BA, 
the subsequent chains are generated using the two previ- 
ous rules, for instance, JF 3 = ABA. Starting with an A, 
we construct the following sequences, A, B,AB, BAB, 
ABBAB, BABABBAB, and so on. Each generation ob- 
tained by iteration of the rules is labeled with an index 
I. Is clear that the number of letters in each generation I 
is given by the Fibonacci numbers F(l) of generation I, 
which satisfy: F(l) = F(l - 1) + F(l - 2) with the initial 
conditions: F(0) = 1,F(1) = l.A Fibonacci potential is 
obtained with two values Va,Vb which follow a FS. 



Our first task in order to compare the FC and the 
Harper potential, consists in finding an analytical expres- 
sion for the Fibonacci potential. This can be done in the 
following way. By using the cut and projection method, 
it is very easy to prove that a one dimensional chain with 
atoms at positions y n determined by a FS are given by 
[24], 

Vn = [n(f>\ 

where the function [_x\ denotes the greatest integer lower 
than x,(f> is a parameter which turns out to be the inverse 
of the golden mean t" 1 = (\/5 — l)/2. The separation 
between atoms is y n +i — Vn- In the Fibonacci sequence, 
the separation takes two values, I a and Is- Thus, we can 
make a Fibonacci potential with the sequence of spacings 
to get, 

V(n) = V B + V A ([(n+l)^- [ivt>\), (3) 

Using the identity x = [x\ +{x}, where {x} is the decimal 
part of x, we obtain that V(n) can be written as, 

V(n) = (V) + 6V({n<j>}-{(n+l)<f>}), (4) 

where (V) = Va4> + Vb{1 — <t>) is an average potential that 
shifts the zero of the energies, and SV is the strength of 
the quasiperiodicity, measured by the difference between 
site-energies 5V = Va — Vb- In what follows, without any 
loss of generality we set Va and Vb in such a way that 
(V) = 0. Now, if the integer variable n is replaced by a 
continuous one, say again x, the resulting potential V(x) 
is just a square wave as shown in figure 1. The potential 
stays at -5V(f> for an interval of length 1 — (j>, and then it 
jumps to 5V(1 — (f>) for a length 4>. 

The decimal part function {x} has period 1, and can 
be developed as a Fourier series, 

1 1 °° 1 
{ x 4>} = 2 ~ ~ ~ sill ( 27r( f>sx), 

s=l 

It follows that, 

oo 

V(n) = V + 25V^V(s) cos (ns<t>(2n + 1)) , (5) 

8=1 

where V(s) is the s harmonic of the Fourier series, V(s) — 
sin(7rs0)/7rs. The first terms of this series are shown in 
figure 1, and in fact, we are approximating a square wave 
by a sum of cosines. Notice the very slow convergence of 
the series due to the 1/s factor of each harmonic. This 
potential can be further reduced if a proper phase \ is 
used in Eq. (4), in such a way that the terms {x<j)} are 
replaced by {x<p + \}- This phase is only a horizontal 
translation of the potential. For \ = —<j)/2, the Fibonacci 
potential is simply written as, 

oo 

V(n) = 26VY,V(s) cos (27rs^n) , (6) 
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FIG. 1: The Fibonacci potential can be obtained by evaluat- 
ing a square wave of period r at integers values. The approx- 
imation with one, two, three and four harmonics, obtained 
from Eq. (6), are shown in the figure. 



It is also worthwhile mentioning that for s — l<f> 1 
(where I is an integer), V(s) — 0. This condition can 
only be hold when 4> is a rational. When this is not the 
case, V(s) rs for integers s an I such that <fi ss l/s, and 
thus l/s is a rational approximant of (f>. Such rationals 
are obtained from the continuous fraction development of 
the number (f>. For the FS, s ~ It from where it follows 
that I and s are successive Fibonacci numbers. Such 
condition can be interpreted as the existence of a rational 
approximant structure when (f> is set to <p = V s - The 
most important Fourier components in Eq. (6), are those 
harmonics for which s i=s (r + l/2)<ft~ 1 for an integer r. 
Using the decomposition in integer and decimal parts, 
this happens whenever {(?' + l/2)(j)~ 1 } is nearly zero or 
one. 

Comparing eq.(2) with eq.(5) and eq.(6), we observe 
that in fact, the Fibonacci potential is just a superposi- 
tion of Harper potentials with the right amplitude and 
phase. In some sense, this is a similar situation to an ap- 
plied effective modulated magnetic field [25]. This leads 
to many questions. The first is how the transition from 
Harper is done. To answer this, we will cut the sum in 
Eq. (6) at a finite number of harmonics, denoted by S. 
In figure 2 a) we plot the energy spectrum of the Harper 
equation in the case A = 1. This spectrum, as well as 
the others discussed in this article, were obtained by us- 
ing the transfer matrix formalism [26]. Figure 1 a) is 
the well known Hosftadter butterfly [3] . Notice that here 
v = n4>, and thus our figure is not exactly the same as in 
Hosftadter work, since therein, a sweep for all values of 
v between and 27r is made, and thus the spectrum has 
more points. 

If (/) is taken as a parameter in eq.(6), then the spec- 
trum of Fig. 2a) is also the same as for the Fibonacci 
potential but with only one harmonic, s = 1. The pa- 
rameter used is SV — tt/ sin(7rr~ 1 ) « 3.3706, chosen to 
correspond to A = 1 in the Harper equation . In Figure 
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FIG. 2: The energy spectrum as a function of the parameter 
(j> for SV — n/ sin(7TT~ 1 ), corresponding to A = 1 in the pure 
Harper equation , using a) one harmonic (Hofstadter butter- 
fly), b) two harmonics, c) three harmonics, and d) the full 
Fibonacci potential ("Fibonacci butterfly"). 



2b) and 2c) we show the effects of adding harmonics s = 2 
and s — 3 in the development of the potential, and fig- 
ure 2d) presents the result for the exact FC. There is an 
important change between the pure Harper case and the 
second harmonic case. Mainly the left part of the Hof- 
stadter butterfly is washed out. Also, is clear that with 
only three harmonics, the structure is already very sim- 
ilar to the pure Fibonacci case. This spectacular wash 
out due to the second harmonic has its origins in the 
difference of lengths in the steps of the square wave Fi- 
bonacci potential. For the first harmonic, this change is 
not observed. But by looking at Fig. 1, one can observe 
that with two harmonics there is an asymmetry in the 
upper part of the series. This effect is more notorious 
when SV > 1, and can be translated in an almost split 
band limit when SV — > oo, around self-energies Va and 
Vb- In fact, this is one of the main differences between 
the Harper and Fibonacci potentials. 

In Figure 3, a similar set of figures presents what hap- 
pens when SV = 1. This corresponds to A = sin(7rr _1 ) sa 
0.2967 in the pure Harper case, and the spectrum is con- 
tinuous. Figure 3b) and 3c) are the cases with two and 
three harmonics. Finally Fig. 3d) shows the case of the 
pure Fibonacci chain, which is a beautiful fractal. Since 
it is known that the FC presents a singular continuous 
spectrum [18], it is open the question for which harmonic 
there is the transition from one type of spectrum to the 
other. 

In all of the previous cases, <j) was studied between 
and 1 since using the following identity, 

cos(27ra;) = cos(27r([a;J + {%})) — cos(27r{a;}), 

we have that cos(2ir(f>sn + k) — cos(2ir {<f>}sn + k). Thus, 
the problem has periodicity 1 in <p, since it depends on <p 
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FIG. 3: The energy spectrum as a function of the parameter 
<f> for SV — 1 using a) one harmonic (corresponding to the 
Harper model at A = 0.2967), b) two harmonics, c) three 
harmonics, and d) full Fibonacci potential. 



only through {</>}. For rational cj) of the form P/Q, with 
P and Q integers, this means P < Q. 



III. DISCUSSION 



In this section we will discuss the main features that 
arises from the previous figures. The spectral properties 
of the pure Hosftadter butterfly has been discussed by 
many others using diverse techniques [3] [4] [27], but in 
order to understand the transition between butterflies, 
here we will explain the main features using the structure 
in reciprocal space of the potential. Let us first study 
the tight-binding equation in the usual approach. We 
propose that the wave- function can be written as [4], 



4>n 



oo 

1 E 

m=— o 



j im(2-K 4>n-\-v) 



where d m is the m Fourier component of the wave func- 
tion. If we introduce this solution into Eq. (1), when 5* 
harmonics are present in the potential, the following set 
of equations are obtained, 

(E* - 2SV* cos(2iT(f>m + k)) d m = d m+1 + c/ m -i+ (7) 

s 

^V*(s) (d m _ s + d m+s ) , 

s=2 

(8) 

where the parameters are defined as E* = 2E/(V(1)6V), 
SV* = 2/{V(l)5V), and, 



V*(s) 



V(s) 
V(l) 



1 sm(nsq 



The other parameter is k = v for the pure Harper equa- 
tion, and k = v = nip for the general case. If only the 
first harmonic of the FC is used, the equation shows 
that a Harper equation in the reciprocal space is also 
a Harper equation with a renormalized set of parameters 
[27]. In the case of Fibonacci, Eq. (7) proves that the 
Schroedinger equation in reciprocal space has a different 
form, since each site interacts with infinite many others. 
The interaction between Fourier components decreases as 
1/s. It is a long range interaction that can be thought as 
a modulating field. Observe that in the transformation 
between the Hofstadter and the Fibonacci butterfly, the 
limitation to <S* harmonics in the potential, is equivalent 
in Fourier space to a cut-off of range S of the interaction. 

When <p is a rational number, say P/Q, the potential 
with any number of harmonics has periodicity Q, as has 
been discussed in the previous section. The correspond- 
ing wave functions are given by, 



e lkn J2 d n 



i) 



(9) 



where the -n/Q < k < ir/Q. The only difference be- 
tween the Harper and the Fibonacci case, is in the values 
of the Q coefficients d m . Since the set d m can be ob- 
tained in the Harper case using quantum groups [8] . It is 
worthwhile mentioning that previous efforts in building a 
perturbation theory to tackle the proprieties of quasiperi- 
odic Hamiltonians have failed due to the small divisor 
problem, that shows up in perturbation expansions. The 
present approach does not have such problem, since a 
quasiperiodic solution can be built from a perturbation 
of a solution that is already quasiperiodic. 

Also, we can show that an effective potential can be 
written in the reciprocal space for periodic approximants 
of the Fibonacci chain. According to Eq. (9), the solu- 
tion must have periodicity Q. with I an integer. Such 
result can also be obtained from V(n) when <f> is the ra- 
tional P/Q- In this case, the factor cos (2irs<f>n) in Eq.(6) 
is repeated for the harmonics s' that has the form s + lQ 
where I is a positive integer. By grouping all the harmon- 
ics module Q, the potential is written as, 



_ Q „ / p 

V(n) = V + 2SV^y r P/Q (s) cos I 27rs-n 



(10) 



where V p/q(s) is an effective potential, 

Osin(7rsP/Q) A^(-l)^ 



Vp,q{s) 



P 



E 

\l=0 



s + lQ 



(11) 



The equation in reciprocal space is reduced as follows, 
(E' - 26V' cos(27T(/)m + fc)) d m = (12) 



s=l 



(13) 
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which shows that for a periodic approximant, the range 
of the interaction in reciprocal space is Q. 

Concerning the spectrum, the band edges are obtained 
from Eq. (7) when k = and k = n/Q. Instead of 
following this path, we will look at how the structure 
of the potential in the lattice reciprocal space determines 
the spectrum for SV < 1. A much more physical insight 
is obtained in this way. This approach is different from 
the one realized in others works [4] [27] , since usually the 
potential is projected in the base e m ( 27 ' ( P n + v ) . Here we 
will project into the reciprocal "vectors" of the lattice (G) 
The main idea is that for a one dimensional crystal, it is 
known that each reciprocal "vector" G with component 
of the potential (V(G)), opens a gap of size, 



V(G) 



at reciprocal vectors q — G/2, G, 3G/2, ... It is possible 
to follow the opening of the gaps by the effect of V{G). 
The reciprocal components are, 



V(G) 



1 



N-l 



£V(n)< 



-iGn 



n=0 



where G can be chosen among the wave vectors q — 
2nt/N, with t = 0, ...,N — 1, in a lattice with N sites, 
were periodic boundary conditions are used. The first 
Brioullin zone is the interval — it < q < n, although to 
simplify the algebra we take for the moment, q between 
and 2tt. Consider first the Fourier component of the 
Harper potential (Vh{G)) for a given parameter 



2A 



N-l 



-i2irtn/N 



(14) 

If ((f) — t/N) is an integer, the problem is almost solved, 
because V(t) = X5((f> — t/N). This happens whenever 
</> = P/Q and, 

P m 

Q = N' 

so when N is chosen as a multiple of Q, as for example 
N = IQ with I = 0, 1, 2, then the harmonic m = P is 
the only one that has a contribution, and the problem is 
solved as a simple cosine potential. 

This solution is very simple compared with other com- 
plex approaches, since what most of the people do is fix 
N and then a sweep of the (f> is made. In other words, 
the Hosftadter butterfly is built for a fixed N. However, 
the present approach shows that if we fix a rational <fi 
and move N for each <fi until N = IQ, the potential is 
much more tractable. In the limit of big I, N ^> Q and 
there are so many vectors k that the condition is very 
easily achieved; in other words, there is a continuum of q 
so e l27 "^ n has a " possible periodicity" of the system. The 
Fibonacci case is a little bit more difficult, because the 



potential is a sum of cosines, but from Eq. (11) is clear 
that the components in Fourier space, that we denote by 
V F (G) , are given by V F (G) = V P/Q (t). 

The problems arise when <j> is an irrational or N is 
not a multiple Q, since the previous trick is not valid. 
Eventually, the sum in Eq. (14) can be made when 4> — 
t/N is not an integer to give, 



V H (G) = V H (t) 



A / 1 - e t2 ^ N 



^ _ e -i2Tr<pN \ 



s/N V 1 ~ e^W-i/W) 1 _ e -i2w(<t>+t/N) 

The corresponding norm is, 

2 A 2 A 



v H (t) 



N B ' 



(15) 



where, 



A = (1 - cos 2iT(j)N) 2 + [cos2tt^ - cos27r^(iV - l)] 2 



2 cos(2TT<f>t/N) [1 - cos 27T0/7V] [ CO s 2ircj> - cos 2-n$(N - 1)] , 



and, 



B = [1 — cos (2tt (4> - t/N))] [1 - cos 2tt (0 + t/N)] . 

Notice how V H (t) -» \6(<j) - t/N) as <p -> t/N. The gaps 
widths depend upon these components, and the maxi- 
mum of Vh (G) occurs when G w ±4>. Now let us propose 
the solution, 



iqn 



for Eq. (1). The resulting equation in reciprocal space 
is, 

(E-2cosq) c q = J2v H (G)c q ^ G 

G 

Since the main contribution to Vh(G) comes from G w 
±(f>, a mixing of wave vectors G « (j> and G w —<p occurs 
at q = ±0/2. Using perturbation theory, this means that 
for A <C 1, the main gaps are open around, 

E w ±Acos(27r(0/2)) = ±A cos(7r0). 

Figure 4 compares this prediction with the Hofstadter 
butterfly, showing an excellent agreement. Other band 
gaps are obtained at q = ±r0/2 with r € N. The general 
pth-ordcr perturbation term is of the form, 

V q = V(q - q\)V(q 1 - q 2 )...V(q N -i - qw) 

where q t labels the vector 2-wt/N. For small A, gaps will 
be open at, 

E w A r cos(7T0r), 

These cosine branches are also plotted in Figure 4, com- 
pared with the Hofstadter butterfly, showing that the 
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This gap is clearly in the horizontal lines at <fi = 1/2 
in Fig. 3 d). The process of folding in reciprocal space 
can be repeated in a similar way for other rationals like 
4> = 2/3 and (p = 1/3. In this case, the periodicity is 3, and 
the folding around the first Brioullin zone limit occurs at 
k = 7r/3. Three bands are produced in this case, and the 
gaps are centered at ±Acos(7r2/3) and ±Acos(7r/3). The 
process is repeated for other approximants.. 



IV. LOCALIZATION PROPERTIES USING 
CORRELATORS 



FIG. 4: General structure of all the obtained spectra. The 
Hofstadter butterfly is compared with the cosine branches 
given by E — Acos(<^m) for m = 1, 2 and 3 for A = 0.2967. 



basic structure of the spectrum is determined by these 
branches. When A is near one, around each gap there 
are many wave vectors that mix together, so the present 
approximation breaks out. 

A detailed observation of the Fibonacci spectrum, 
shows that the main effect upon the spectrum is an asym- 
metry of the cosine branches, and the different gap sizes 
around these branches, due to the mixing of different 
modulating harmonics. This is very clear in Figures 3 
b) and 3 c). The same analysis performed to the Harper 
equation is valid for these cases, specially for low S. . 
The equivalent rough approximation for the FC is, 



E : 



sin(7rs</ 

ITS 



C0S(7T<J 



What is behind the process of gap opening, is a self- 
similar folding of the bands in the reciprocal space, as 
revealed by following the sequence of approximants for <p 
in all the obtained spectra. For <p — 0, the spectrum is 
continuous and goes from -25V to 25V. The next most 
simple spectrum corresponds to <f> — 1/2, which gives 
two bands. If the lattice has period 2, then the first Bril- 
louin zone of the zero approximant (4> = 0), is folded 
around half the original Brioullin zone limit, k = tt/2. 



V(G] 



is open 



As usual, an energy gap of size Aq ~ 2 

in the zone boundary due to a mixing of waves with re- 
ciprocal vectors G/2 and -G/2. In this case, G =tt/2 
and V'(G) = 2 A. For the Harper equation we get, 



A 



G 



4A, 



In the Harper equation, there is symmetry around E = 0, 
thus the bands limits are E = ±2A. For the pure FC first 
approximant {<p = 1/2), V{G) produces a different gap, 
and the central gap limits are, 



As an example of the utility of having an expansion 
of the Fibonacci potential, we show how to obtain cer- 
tain localization properties from the potential correla- 
tors. This can serve to understand how the addition of 
harmonics leads to different localization properties. For 
A = 1 the eigenstates of the Harper equation are critical, 
while they arc localized for A > 1 and A < 1 . This comes 
from Eq.(7). Therein, if the square of the components is 
finite. 



E 



\d„ 



< oo 



and the corresponding wave-function is non-localized. In 
the pure Harper equation, when A — > oo, the solutions 
in real space are localized wave-functions, since it cor- 
responds to the dual of a solution A — > in reciprocal 
space, which are known to be extended. 

For the Fibonacci chain, it is known that for all values 
of A the eigenstates are critical [26]. An interesting ques- 
tion is how many harmonics are needed to produce this 
transition. Although this requires a complete study, let 
us show how this problem can be tackled. To measure 
localization in one dimension, the Lyapunov exponents 
are used {L~ r ). They are the inverse of the localiza- 
tion length (L) for an exponential localized state. It has 
been shown that for one dimensional Hamiltonians, the 
localization properties depend upon the pair correlation 
of the potential [28] when V(n) <C 1. Notice that this 
is the most interesting limit for quasipcriodic potentials, 
since for V[n) 3> 1 there are already many methods to 
treat the spectrum and the eigenfunctions [29] . For a one 
dimensional Hamiltonian [28], 



L = 



8sin 2 0) : 



2^£(fc) cos(2/ifc) 



k=i 



±5V 



Here, the function <p(fi) is given by the Fourier series 
with the coefficients £(fc) which is the correlator of the 
site potential V(n). The correlator £(fc) is defined as, 

(V(n)V(n + k)) = efc(fc), (V{n)V(n)) = el 

For the Harper equation, £(fc) = cos(27r0fc), which 
gives ip(n) = and thus all states are non exponentially 



7 



localized. In the case of a FC, by using Eq. (6), the 
correlator can be written as, 



C(fc) = 45V£ 



sin 2 (7TS() 



■ C0S(27T S(f)k) 



which also gives (f(fJ.) = for all S. Thus, we have proved 
that all states arc non exponentially localized for any 
number of harmonics. The states are extended or criti- 
cal. An interesting question that remains to be answered, 
is how the transition from extended to critical states is 
achieved as the number of harmonics is increased from 
Harper to Fibonacci. To solve this question, an expres- 
sion for the scaling exponents in terms of the correlators 
is needed [30] . 



potentials. As a consequence, one can follow the evolu- 
tion of the spectral types as a function of the number of 
harmonics. In particular, a butterfly similar to the Hofs- 
tadter case is found for the Fibonacci. The corresponding 
spectrum is a fractal object. The Fourier components of 
the potential provides a simple explanation for the main 
features of the spectra. We also show how the devel- 
opment of the potential leads to interesting questions, 
as for example, at which harmonic the spectral type is 
changed. This is equivalent to ask for which harmonic 
the eigenf unctions of the Fibonacci case become criti- 
cal. The present approach also leads to the possibility 
of building analytical solutions for a FC. We hope that 
other researchers will answer some of these interesting 
questions in the nearby future. 



CONCLUSIONS 



In the present article, we have shown that the Fi- 
bonacci potential can be written as a sum of Harper 
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